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We derive a spacetime formulation of quantum general relativity from (hamiltonian) loop quan- 
tum gravity. In particular, we study the quantum propagator that evolves the 3-geometry in proper 
time. We show that the perturbation expansion of this operator is finite and computable order by 
order. By giving a graphical representation a la Feynman of this expansion, we find that the theory 
£j>^ | can be expressed as a sum over topologically inequivalent (branched, colored) 2d surfaces in 4d. The 

contribution of one surface to the sum is given by the product of one factor per branching point of 
the surface. Therefore branching points play the role of elementary vertices of the theory. Their 
value is determined by the matrix elements of the hamiltonian constraint, which are known. 

The formulation we obtain can be viewed as a continuum version of Reisenberger's simplicial 
quantum gravity. Also, it has the same structure as the Ooguri-Crane-Yetter 4d topological field 
theory, with a few key differences that illuminate the relation between quantum gravity and TQFT. 
■ Finally, we suggests that certain new terms should be added to the hamiltonian constraint in order 

to implement a "crossing" symmetry related to 4d diffeomorphism invariance. 
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Q . I. INTRODUCTION 



An old dream in quantum gravity [jD is to define a manifestly spacetime-covariant Feynman-style "sum over trajec- 
tories" O , sufficiently well defined to yield finite results order by order in some expansion. The hamiltonian theory has 
obtained encouraging successes in recent years, but it suffers for the well-known lack of transparency of the frozen time 
^ formalism, for the difficulty of writing physical observables and for operator ordering ambiguities. These problems are 
I \ related to the lack of manifest 4d covariance. Here, we derive a covariant spacetime formalism from the hamiltonian 
£~ I ■ theory. This is of course the path followed by Feynman to introduce his sum over trajectories in the first place [[|. 
7v*| What we obtain is surprising: we obtain a formulation of quantum gravity as a sum over surfaces in spacetime. The 
J> . surfaces capture the gravitational degrees of freedom. The formulation is "topological" in the sense that one must 
sum over topologically inequivalent surfaces only, and the contribution of each surface depends on its topology onlyf). 
This contribution is given by the product of elementary "vertices" , namely points where the surface branches. The 
sum turns out to be finite and explicitly computable order by order. The main result of this paper is the construction 
of this finite "sum over surfaces" formulation of quantum gravity. 

Let us sketch here the lines of the construction. Given gravitational data on a spacelikc hyper-surface £j, the 
three-geometry on a surface £ / at a proper time T in the future of £ j (as measured along geodesies initially at rest on 
£,•), is uniquely determined in the classical theory. It is then natural to study the corresponding evolution operator 
U(T), that propagates states from £j to £/ in the quantum theory. This operator, first considered by Teitelboim 
PL codes the dynamics of the quantum gravitational field, and is analogous to the Feynman-Nambu proper time 
propagator || for a relativistic particle. Here, we construct the operator U(T) in quantum GR, and we expand it in 
powers of T. We obtain a remarkable result: the expansion is finite order by order. This is our first result. 

Next, we construct a graphical representation of the expansion. This is obtained by observing that topologically 
inequivalent colored 2d surfaces a in spacetime provide a natural bookkeeping device for the terms of the expansion. 
We obtain an expression for U(T) as a sum of terms labeled by surfaces a bounded by initial and final states. A 
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1 More precisely: on the diff-invariant properties of the surface. 
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surface a consists of simple components (2-manifolds), or "faces", that carry a positive integer or color. Faces meet on 
Id "edges" , colored as well. Edges, in turn, meet at branching points, denoted "vertices" . The weight of each surface 
in the sum is a product of factors associated to its vertices. The value of a vertex is determined by the hamiltonian 
constraint, and is given by a simple function (involving Wigner 3n-j symbols) of the colors of the adjacent faces and 
edges. This "sum over surfaces" version of the dynamics of quantum general relativity is our second result. 

The construction allows us to consider transformation properties of the hamiltonian constraint under 4d diffeomor- 
phisms (diff) in a manifestly covariant way. This analysis suggests the addition of certain new terms to the constraint, 
corresponding to an alternative operator ordering, which implements a "crossing" symmetry at the vertices, with a 
nice geometrical appeal. Thus, 4d diff invariance may be a key for reducing the present ambiguity in the operator 
ordering. Furthermore, the new terms seem to prevent some potential problems with locality pointed out by Smolin 
. The introduction of these new terms in the quantum hamiltonian is our third result. 
The idea that one could express the dynamics of quantum gravity in terms of a sum over surfaces has been advocated 
in the past, particularly by Baez |?|] and Reisenberger ||. On the lattice, a sum over surfaces was recently developed 
by Reisenberger and the lattice construction has guided us for the continuum case studied here. It is important 
to emphasize, however, that the present construction is entirely derived from the canonical quantum theory in the 
continuum. 

The sum over surfaces we obtain has striking similarities with topological quantum field theory (TQFT). More 
precisely, it has the same kinematic as the Ooguri-Crane-Yetter model a 4d TQFT which extends the Ponzano- 

Regge-Turaev-Viro 3d TQFT [jl2] [uj to four dimensions. Essentially, the difference is given just by the weight of the 
vertices. In Appendix B we discuss similarities and differences between the two theories. The discussion, we believe, 
sheds much light over the tantalizing issue |l5| of the relation between finite-number-of-degrees-of-freedom TQFT and 
quantum gravity. In particular, we argue that a diff-invariant quantum field theory with an infinite number of local 
(but non-localized) degrees of freedom -such as quantum general relativity- can be obtained by having a sum over 
arbitrarily fine triangulations, instead of triangulation independence, as in combinatorial TQFTs. 

On the other hand, the sum over surfaces we obtain can be viewed as (a first step towards) a concrete implementation 
of Hawking's sum over 4-geometries 0. In fact, the surfaces over which we sum have an immediate interpretation 
as "quantum" 4-geometries, as we will illustrate. This fact should make the the general techniques of covariant 
generalized quantum mechanics fllf | available to quantum gravity, potentially simplifying the difficulties with physical 
observables of the hamiltonian formalism. 

The basis of our construction is loop quantum gravity []l7|-[l9||. The finiteness of the sum-over-surfaces and the 
picture of a "discrete 4-geometry" that emerges from this work are related to the fact that geometrical operators have 
discrete spectra. The discreteness of the spectra of area and volume -and the "quantized" structure of space that 
these spectra suggests- is a central result in loop quantum gravity, first obtained by Rovelli and Smolin in pp| |, and 
later confirmed and clarified by a number of authors [^l] . The main ingredient of our construction is the quantum 
hamiltonian constraint |l7]|32]]. In particular, Thiemann's version of the hamiltonian constraint |33| and some variants 
of it play an essential role here. Matrix elements of this operator have been computed explicitly in |Q, using the 
methods developed in [^jf! 

In section II, we summarize the basics of nonperturbative loop quantum gravity. In section III we define the proper 
time propagator and its expansion. In section IV we show that the proper time propagator can be expressed as a 
sum over surfaces. In section V we discuss crossing symmetry and the new terms of the hamiltonian constraint. 
In section VI we summarize and comment our results. Appendix A is a brief glossary of some geometrical term 
employed. Appendix B contains the comparison with TQFT. In Appendix C we give an example of 3d diff-invariant 
scalar product. A preliminary version of this work has appeared in fl3l| . 



2 Note added to the revised version: An extremely interesting sequence of developments has appeared after the first version of 
this work. Smolin and MarkopoulouJ 2& have attempted to incorporate the causal structure of the lorentzian theory directly 
into the surfaces' kinematic. Baez |27 1 has explored the general structure of the theories of the kind introduced here, namely 
defined by a sum over surfaces with weights at the vertices. Markopolou jii) and Baez J2s| have illuminated the nature of the 
vertices by consider evolution in the context of triangulated spacetimes. In a triagulated manifold "time evolution" can be 
decomposed into a finite number of elementary moves, corresponding to the possible ways of adding a (spacetime) (3+l)-simplex 
to a (space) 3-d boundary. When viewed in the dual/ of the triangulation, these elementary moves correspond to the (crossing 
symmetric) vertices generated by the hamiltonian constraint operator and studied in this paper. We find this convergence 
fascinating. Finally, Crane ]3^] and Baez J2i| have explored two possible ways of obtaining quantum GR by modifying the 
Ooguri-Crane-Yetter state sum. (On this, see also |^|.) We think that deriving the sum over surfaces formulation of quantum 
GR directly from the lagrangian would be a major step ahead and would shed light on the problem of the integration over the 
lapse. 
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II. CANONICAL LOOP QUANTUM GRAVITY 



A. Kinematic 

We start with nonperturbative canonical quantum gravity in the loop representation . The Hilbert space Ti of 
the theory is spanned by the basis | S) , where S is a spin network []35| . A spin network is a colored graph Y embedded 
in a (fixed) three dimensional compact manifold For a fixed choice of a tree- like expansion at the nodes, these 
states are orthonormal p3| , p^6|j37| (see |23| for details on their normalization). 

(S'\S) = 6ss>; (1) 

the matrix elements of the change of basis between different nodes' expansion can be derived from equation (|C7|). An 
equivalent construction of this Hilbert space can be obtained in terms of functions over (generalized) connections jl8) . 
For details on the equivalence between the two formalism, see |3j| . 

The dynamics of quantum general relativity is governed by two operators in Ti.: the Diff constraint operator C[N] 
and the hamiltonian constraint Cl [N] . Let us examine them. 



B. The Diff constraint and its solutions 



For every diffeomorphism / : X — > E (in Diffo, the component of the diffeomorphism group connected to the 
identity), let D[f] be the operator in Ti giving the natural action / : S i — ► / ■ S of the diffeomorphism on the spin 
network states. Namely 

D[f]\S) = \f-S). (2) 

For every vector field N on £ that generates a one parameter family f t of diffeomorphisms (by dft/dt = N, 
/o=identity)Q the Diff constraint is defined by^| 

; (3) 

t=o 

and corresponds to the classical diffeomorphism constraint smeared with Shift function N. The space Tidiff of the 
solutions of the diffeomorphism constraints is defined as H-diff = Diff ' ^ * s s P anne d by a basis s), where s is 
an s-knot, namely an equivalence class of spin networks under diffeomorphisms, which define the linear structure of 
'Hdiff- One can define the scalar product in T~Ldiff by an integration [ ^Opl] ] over Diffo. li S G s and S' G s', 

(s\s')=Af[ [df}{f-S\S'). (4) 

J Diffo 

Af is a normalization factor. Equation (Ji|) is meaningful because the integrand vanishes over most of the integration 
space (because two spin network states are orthogonal unless they have the same graph) and is constant on a discrete 
number of regions whose volume is normalized to one by J\f . Thus we have 



C[N] = -i- D[f t ] 



3 We recall the definition of coloring of a spin network j2^] . Each node of the graph with valence higher than three (more than 
three adjacent links) is arbitrarily expanded in a tree-like trivalent subgraph. The internal links of the subgraph are denoted 
virtual links. The coloring of the graph is an assignment of a positive integer to each real or virtual link - in such a way that 
at every trivalent node the sum of the three colors is even and none of the colors is larger than the sum of the other two. The 
set of the colorings of the virtual links of a node is also called coloring of the node. A coloring can be thought as an assignment 
of an irreducible 517 (2) representation to each link and of an invariant coupling tensor to each node. 

4 We put an arrow over vectors (N); but not over spatial coordinates (x) or diffeomorphisms (/). 

5 Rigorously speaking, C[N] is not well defined on Ti. This is due to funny (kinematical) inner product (|l|), in terms of which 
the action of the diffeomorphism group is not strongly continuous. This fact does not disturb the construction of the theory, 
because the only role played by C[N] is to implement invariance under the finite transformations it generates. These are well 
defined |3q|. Here, it is useful to consider C[N] as well, because it plays a role in the formal manipulations below. 
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(8\s f ) =Y,(pSi\S'), 
p 



(5) 



where the sum is over the (discrete) automorphisms p that send the graph and the links' coloring into themselves. 
See Appendix C for an example, and reference for a rigorous construction. It is useful to view an s-knot state 

as a group integral of a spin network state: 

H=W W](f-s\, (6) 

J Dtffo 

where S G s. 



C. The hamiltonian constraint 

The Hamiltonian constraint that we consider is the density-weight 1 hamiltonian density, smeared with a density- 
weight Lapse function N. The Lorentzian hamiltonian constraint Cx[7V] can be written as the sum of two terms: 
Cl[N] = C[N] + V[N] [Q, where Cl[N] is the Euclidean Hamiltonian constraint. For simplicity we deal here only 
with the first term. Thus, we are dealing below with Euclidean quantum gravity only. We expect the methods 
developed here to be extendible to the V[N] term as well, and therefore to Lorentzian GR, using the techniques 
developed by Thiemann ||| . 

The definition of C[N] is plagued by ordering ambiguities [^j33l|43fl . Some of these are fixed by 3d diff invariance 
p2| . In section V we discuss how 4d diff invariance might fix others. Here, we recall Thiemann's version of the 
Hamiltonian constraint, which is the starting point of our construction. First, the non-symmetric operator C ns [N] is 
defined as 

C ns [N]\S) = 2 J2 N ( x i) £ Yl A lJKt AS) D lJKtt , \S). (7) 

i£n(S) (Jif)Ge(i) e=±l,e'=±l 

Here i labels the nodes of S (which form the set n(S)); Xi are the coordinates of the node i; (JK) labels the couples of 
distinct links adjacent to the node i (these form the set e(i)); the operator DijKee/ was introduced in f43]| , it acts on 
the spin network by creating two new trivalent nodes %' and i" on the the two links J and K respectively, connected 
by a link with color 1, and adds e (resp. e') to the color of the link connecting i and i' (resp. i and i"). This is 
illustrated in Figure 1. 




FIG. 1. Action of D iJKee i . r, q and p are the colors of the links I, J and K. 

The precise location of the nodes and the link added is an arbitrary regularization choice. The coefficients Auj tt i (S) 
of Thiemann's operator are well defined and can be computed explicitly |54j . They are functions of the colors of the 
links adjacent to the node i; they are finite and can be expressed as products of linear combinations of n-j symbols 
of 517(2). 

It is important to notice that Thiemann's operator was derived using the infinite-dimensional differential-geometry 
techniques introduced in |l8| . These differential techniques were introduced as a mathematical sistematization of the 
ideas on loop quantization introduced in p7| |. They have shed much light on loop quantum gravity, have provided a 
rigorous mathematical foundation of the theory, and have led, among other results, to Thiemann's operator. However, 
the operator itself is a well-defined algebraic operator on the spin network basis, and the computation of its matrix 
elements is easier using algebra than using infinite-dimensional differential-geometry p3| . The equivalence between 
the purely algebraic formalism ( "loop representation" or "spin network representation" ) and the differential formalism 
( "connection representation" ) is shown in detail in (39| (see also ||^],^| ) . The situation is analogous to the two well- 
known ways of computing the spectrum of the harmonic oscillator: one can use Dirac's algebraic technique, in the 
\n) basis; or, alternatively, one can solve the differential Schrodinger equation, in the coordinate basis. Each basis 
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has its own advantages, but the two formalisms are equivalent, and there is no sense in which one representation is 
more "rigorous" than the other. In particular, spectra of area and volume in quantum gravity can be computed cither 
using algebraic techniques (this is the way they were first computed and their discreteness was discovered in or 
using differential techniques (see |25fl , and references therein). The resulting spectra are, of course, equal j24]. 
We shall use Thiemann's symmetric hamiltonian constraint operator C[iV], defined by symmetrizing C ns [N] 



(S'\C[N]\S) = ±({S'\C ns [N]\S) + (S\C ns \S')). (8) 



Explicitly, we have 



(S'\C[N]\S)= J2 N ( X i) E A iJK ee,(S) (S'\DiJKee,\S) 
iGn(S) (JK)eel 

+ E N ( x *') E A ilJKe , e ,(S)(S>\Dt JKee ,\S) (9) 

i'en(S') (JK)eel 

where 

(S'\D\ JKeel \S) = (S\D iJKee ,\S'). (10) 

An observation important for what follows is that the hamiltonian constraint is "local" , in the following sense. 
Given a spin network S, we may cut it in two parts, by cutting n links, obtaining two spin networks with open 
ends, which we denote S and S. Imagine we have two spin networks Si and Sf that can be cut as Si and Si, and, 
respectively, Sf and Sf. Imagine that Si = Sf, namely the two spin networks differ only in their "tilde" component. 
Then the matrix elements {Sf \C[N] \Si) do not depend on the "hat" components Si and Sf, so we can write 

(St\C[N\\Si) = (S f \C[N]\Si). (11) 

This decomposition will play an important role below. 

We simplify notation by introducing a single discrete index a, (3 ... to replace the discrete set of indices (i, JK, e, e'). 
For every spin network S, a ranges over a finite set [S] of values, with 

4 n (12) 

iGn(S) 

values, where Vi is the valence of the node i. We also indicate by x a the coordinates of the node with index a. Using 
this, we have 

(S'\C[N]\S) = N(x a ) A a (S) (S'\D a \S) + £ N(x p ) A P (S) (S"|Dj|5>. (13) 

ce[S] 0e[S>] 

The hamiltonian constraint transforms covariantly under the diffeomorphisms generated by the diffeo constraint 

'C[N],C[N]]=C[£^N], (14) 

where is the Lie derivative along N . Under a finite diffcomorphism /, we have 

D[f] C[N] D- l [f]=C[N f ], (15) 

where Nf is the transformed Lapse: 

Nf(x) ee N(f(x)). (16) 

The transformation properties of C[N] under 4d diffeomorphisms are less clear. In the canonical formalism these are 
controlled by the commutator of C[N] with itself, which, however, is not fully under control, due to the interplay be- 
tween regularization and 3d diff invariance p^] . (Notice that in 1 33 1 it is shown that the commutator [C„.s[iV], C ns [M]] 
vanishes on diff invariant states; the commutator [C[iV], C[M]] is more tricky.) In Section IV we suggest a way for 
addressing these difficulties. 
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III. PROPER TIME EVOLUTION OPERATOR U(T) 



Consider, as an illustrative example, the Schrodinger equation for a single particle in a potential. If H is the 
Hamiltonian operator, the equation is formally solved by the evolution operator 

U(t) = U(t, 0) = e _< So dt ' H{t ' ] (17) 

where exponentiation, here and below, is time ordered. The matrix elements of this operator between position 
eigenstates define the propagator 

P(x,t-x,t r ) = (x\U(t,t')\x), (18) 
from which the solution of the Schrodinger equation with initial data ip{x' , t') at t' can be obtained by simple integration 

iP(x,t) = J dx'P{x,t;x't')tp{x',t'). (19) 

Under suitable conditions, the propagator can be computed by means of a perturbation expansion in the potential, 
and the expansion has a nice graphical representation. 

For a (free) relativistic particle, we have the option between using the above formalism with the relativistic Hamilto- 
nian (H = \/ p 2 + m 2 ) , or using a manifestly covariant formalism. This was originally done by Feynman by changing 
the description of the dynamics: instead of representing motion by means of the evolution of the three variables x in t, 
we consider a ("fictitious") evolution of the four variables x — (x, t) in the proper time T. This evolution is generated 
by the operator H = p 2 — m 2 = (p ) 2 — p 2 — m 2 . The corresponding proper time evolution operator and proper time 
propagator are 

U(T) = e ~ i So dt ' H{tl) (20) 

and 

P(x,t;x',t'; T) = (x, t\U(T)\x", t'). (21) 

The relation between this proper time propagator and the physical propagator (which is the quantity we compare 
experiments with) is given by 



/>oo 

P(x, t; x', t') = dT P(x, t; x' , t'; T). 
Jo 



(22) 



Or 

P(x,t;x',t') = (x,t\U\x , ,t / ), (23) 

where 

/•OO 

U = dT U{T). (24) 

This can be verified by means of a simple calculation. U is the projector on the physical state space, which codes the 
theory's dynamics. 

Alternatively, one can consider evolution in a fully arbitrary parameter t. Such evolution is generated by H (t) = 
N(t)H, where N(t) is an arbitrary Lapse function. The corresponding evolution operator is 

U N ^ e - l So dtN ^ H , (25) 
which is related to the physical U by the functional integral 

U = j \dN] U N . (26) 

This functional integration can be split into two parts by defining the proper time T in terms of the Lapse as 



G 



T = [ N(t) dt. (27) 
Jo 

Using this, we can first integrate U n over all Lapses N having the same T. 

U(T) = / dN U N (28) 



JT 

where the subscript T indicates that the functional integral must be performed over all TV's satisfying [2?]. Then we 
integrate over T to get the physical quantity U. Thus we have 

U N i — > U(T) i — > U. (29) 

An important observation is that in order to compute the functional integral ^8| we can simply gauge fix N, requiring, 
for instance dN(t)/dt = 0. In this gauge, the integral becomes trivial, and we have U{T) = U(T), which is given in 
|20| . In fact, the functional integration over TV is largely trivial, since Un depends on N only via T. 

We are now going to follow the same path in general relativity. In particular, we will concentrate here on the 
definition and the computation of the proper time evolution operator U(T) and the corresponding proper time 
propagator (its matrix elements) for quantum general relativity. 



A. Definition and meaning of the proper time propagator in general relativity 

In the canonical theory, the ( "unphysical" or coordinate-) evolution of the gravitational field is generated by the 
hamiltonian 



h n #(*) = / Ml X ) C ( X ) + # °(*> x)C a {x)\ = C[N(t)] + C[N{t)] (30) 



a 

(units are fixed here by h — c = IdirG Newton = 1, and we take £ compact). The quantum evolution operator that 
evolves from an initial hyper-surface Sj at t = to a final hyper-surface at t = 1 is 

U N j = e~ l fo dt ff «« (t) . (31) 
We define the proper time evolution operator for quantum gravity as 

U(T) = f [dN, dN} U N # , (32) 

where the subscript {T, *} means that the integral is over all Shifts and Lapses that satisfy 

N(x,t)=N(t) (33) 

i 

dtN(x,t)=T. (34) 

Notice that T is the proper time separation between Ej and £/, defined as the reading on £j of the free falling test 
clocks that started off at rest on . This is because if the Lapse is constant the geodesies that define the proper time 
foliation remains normal to the ADM hypersurfaces. We denote the matrix elements of the operator (|32|) 

P(s f , Si ;T) = (s f \U(T)\ Si ) (35) 

as the proper time propagator [Q. In this paper, we focus on this quantity. We compute it as a power expansion in 
T in the next subsection, and show that it admits a sum over surfaces representation in section IV. n 



Generally, slicing by equal proper time hypersurfaces develops singularities, because the geodesies that define the proper time 
intersect. This causes the canonical evolution to break down: a coordinate system based on the slicing develops coordinate 
singularities where the ADM momentum density diverges. We ignore these difficulties here, but two comments are in order. 
First, the explicit expression for U(T) that we obtain is simple and well defined order by order for any T. We think that 
potential singularities in U (T) should be looked for directly in the quantum formalism. Second, Lewandowski |^ has pointed 
out that the classical evolution of the Ashtekar's variables is better behaved than the ADM variables, because Ashtekar's 
variables can be represented as differential forms, whose components are well behaved at coordinate singularities of the type 
being considered. 
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The construction generalizes to a multifingered proper time. In this case, let E/ is given by t — t(x). The coordinate 
time evolution operator from Ej to T,f at fixed Lapse and Shift is 



jj _ = e ~i f s d 3 x j Q W dt [N{t,x)C(x)+N a {t.x)C a {x)\ ^ 

The multifingered proper time evolution operator is 

U[T]= f [dN, dN} U NI j , (37) 

J[T] 

where the subscript [T] means that the integral is over all Shifts and Lapses that satisfy ( |33| ) and 

r*(9 t _1 (*)) 



dtN(g- 1 (x),t)=T(x) (38) 



o 



where g t is the finite, time-dependent, transformations of spatial coordinates generated by integrating the shift: 

90 (x)=x ) ^M = N(x,t), g = 9l . (39) 

gt and g are functionals of the Shift. As before, T{x) gives the proper time separation between the two hyper-surfaces, 
defined as the proper distance of E/ from Si along a geodesic starting at rest on E$ on x Indeed, if the Lapse is spatially 
constant, the geodesic will be (x) (because in the coordinates (y, t) = (g t (x),t) the Shift vanishes, the Lapse is 
still constant and therefore y — constant is a geodesic normal to all ADM slices). And therefore the geodesic that 
starts off in x reaches E/ at the time t determined by t = t(g^ 1 (x)). Notice that the two notations U[T] and U(T) 
indicate different objects. U[T] is the multifingered proper time propagator, a functional of the function T{x)) while 
U(T) is the proper time evolution operator, which is the value of U[T] for T{x) — constant — T. 

In the rest of this subsection, we discuss the physical meaning of the quantity we have defined and its role in the 
theory. 

First of all, the proper time propagator is the first step toward the computation of the physical evolution operator 
U, in the same sense as the Feynman-Nambu proper time propagator. The operator U is given by functionally 
integrating U N jf over all lapses and shifts: 

U = J [dN] J [dN] U N j. (40) 

This functional integrations corresponds to the implementation of the canonical constraints. As for the relativistic 
particle considered in the previous section, we can split the computation of U from U N ft into two steps 

U N ^ — » U[T] — > U, (41) 

by first computing the propagator U[T] at fixed T(x) and then integrating over T(x). As for the particle, we can 
partially fix the gauge in which we compute t/[T]. In particular, we can choose to integrate over spatially constant 
lapses only. Therefore we have U[T] = U[T], which is given in|37]. And, as for the particle, we can write 



U = J [dT] U[T] (42) 



Thus, U is just the integral over proper time of the multifingered proper time evolution operator j37|). 

We add a general argument that better illustrates why we can fix a gauge for computing J7[T]. This argument is 
formal, but it is interesting because it illuminates the relation between what we are doing and the sum over geometries 
considered by Hawkings 0] . In the metric formulation of canonical GR, U can be written as a sum over four- metrics 
bounded by given initial and final three-geometries. Each such four-metric determines a proper time separation T(x) 
between the initial and final hypersurfaces. Therefore, the integral can be split in two parts, first the integral U(T), 
restricted to 4- metrics with total multifingered elapsed time T{x), then the integration over T{x). In computing U (T), 
we can change integration variables from the four-metric to the ADM variables, namely three-metric, Lapse and Shift. 
The integral contains a high redundancy, corresponding to diffeomorphism gauge-invariance (as the corresponding 
integral for the particle in the previous section did). We can fix part of this redundancy with a condition on the 



8 



Lapse. If we pick an arbitrary Lapse, the condition that the proper time separation of the initial and final slices is 
T{x) becomes a condition on the three-metrics over which we are integrating. However, we can choose a spatially 
constant Lapse, satisfying [33] and [38|. We can always do that, because we may always slice a four geometry with 
equal proper time hypersurfaces, with the result that the corresponding Lapse and Shift satisfy the ( |33"1) and (pSj) . 
Conversely, any history of 3-metrics, together with a Lapse and Shift satisfying (^3|) and ( |38| ) defines a 4-geometry 
with elapsed proper time T{x). In this way, we implement the T(x) condition, without interfering with the integration 
over three-metrics. Thus, U[T] can be computed by fixing a Lapse satisfying (^3|) and (|38|). 

The operator U is a key quantity for the theory. Computing it virtually amounts to solving the quantum constraints, 
including the hamiltonian constraint, which codes all the dynamics of the theory. There are various of ways of looking 
at U. First of all, it is the projector on the the physical state space of the theory. Second, the scalar product (s\U\s') 
defines the physical scalar product of the theory. Therefore 

(S, ^physical = J [dT] (s\ U[T] | s'). (43) 

Finally, we can view matrix elements of U as observable transition amplitudes between quantum states. The details 
of the interpretation of U will be discussed elsewhere, but in all these instances, the role of U is just analogous to its 
counterpart for the single particle. (In this paper we do not attempt to compute U .) 

Does the proper time propagator P(sf,Si; T), have a direct physical interpretation? A simple answer is that 
P(s-f,Si; T) codes the dynamics of the theory, but it has no direct physical meaning: only after integration over 
proper time we obtain a quantity that we can, in principle, campare with experiments. 

This said, we can nevertheless assign a plausible physical interpretation to the proper time propagator P(sj, s^; T), 
with a certain caution. This would be particularly useful for helping intuition. Let us return to the rclativistic 
particle. In that case, Feynman considered the "fictitious" evolution of x and t in the proper time T. Classically, this 
is not incorrect (because the equations of motion of the 4+1 dimensional theory give the correct physical 3+1 motion) 
provided that one remembers that the degrees of freedom are 3 and not 4. Quantum mechanically, in the fictitious 
theory we are quantizing one variable too much. Taken literally, the particle proper time propagator describes the 
3 degrees of freedom of the particle position, plus the degree of freedom of an extra quantum variable sitting on the 
particle, growing with proper time, and not affecting the particle's motion^]. With such a little clock on the particle 
(say the particle is an oscillating molecule) , we could make experiments we could compare the proper time propagator 
with. 

In general relativity, such a "fictitious" evolution with extra degrees of freedom is provided by the so-called "local 
interpretation" of the theory (see [[l6| for a detailed discussion). In this interpretation, the coordinates are interpreted 
as labels of reference-system physical objects (RS-objects). It follows that local quantities are physical observables, 
and that the lack of determinism of the Einstein equations can be interpreted as a consequence of the fact that 
the dynamical equations of the RS-objects are neglected. Under this interpretation, GR is approximate (because 
we disregard the RS-objects energy momentum) and incomplete (because we disregard the RS-objects dynamical 
equations). The incompleteness leads to the apparent physical indeterminism. If we adopt this view, then we can say 
that s, s' and T are observable, because E^ can be physically specified by the RS-objects, and we can use RS-clocks 
to find out where is In doing so, we take approximations that might be ungranted (on the quantum behavior of 
the RS-objects). Concretely, one may consider a definite model, for instance the "dust" model introduced in ||(| and 
studied in [47|. In such a model, a E^ to Ey propagator (where E^ is defined dust variables) is an observable quantity. 
In a suitable limit in which the dust physical effects are disregarded, such propagator might be approximated by the 
pure gravity proper time propagator P(s/, T). 

With all these caveats, one can intuitively think of P(s/,Si; T) as the quantum amplitude that the quantum 
gravitational field be in the state s/, T seconds after being in the state Si. 

B. Expansion of the proper time propagator 

We begin with an observation. Let us write U N ^ as a limit of products of small time propagators. Writing e = 1/K 
and tk — fee, for integers K and fc = 1...K, we have 



7 More precisely, the extra variable is not the evolution of this proper-time clock- variable T in a Lorentz time a; , but rather 
the evolution of x° in T. 
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U - = lim e~ le - f/ «^ (tK) e - ieH N,N(^) e ~ iiH N,N^) 

N ' N K^oo ' 

= lim e -^C[N(t K )] e -ieC[N(t K )] _ _ _ e ^ C [N(t 2 )] e ~ieC[N(t 2 )\ e -ieC[N{tJ] e -ieC[N{t^)\ 
K^oo 

= lim e- ieC ^ tK ^D[f K ]... e ~ ieC ^ D[f 2 ] e~ le c[Nt > ] D[h] 

= lim^e-^"^ D[f K ] . . . ({D[f 3 ]D\f 1 ])- 1 e- i ' ^^D\f2]D[f 1 ]) (p-^e-^^ D[fi 
= lim D 

K~*oo 



[g] (D-i[g K }e-^ N ^D[g K }).., (D-'[g 2 }e'^ N ^D[g 2 }) (D-^e-W^Dfc 



lim D[g) e ~ ieC[A, « (tK)1 . . . e - ieC{N ^ t%)] e- ieC[N ^ tl)] 
= D[g] U Nft , . (44) 



Here is defined by N(x,t) — N(g t (x),t), namely it is the lapse in the coordinates obtained by integrating the 
shift, ft is the "small" diffeomorphism generated by the Shift between the slices t — 1 and t; while gt is the finite 



diffeomorphism generated by the Shift between the slices t — and t. The first equality in (44) is just one of the 
definitions of the time ordered exponential. The second is based on the fact that for sufficiently small time interval 
1/K (sufficiently high K) one can disregard the commutator term in disentangling the exponent (this term is quadratic 
in 1 / K) . The third equality is simply a rewriting of the exponent of an infinitesimal diffeomorphism as a finite (but 
"small") diffeomorphism. The fourth equality is simply the insertion of terms like (Z3[/i]Z? _1 [fi]) in suitable places. 
The fifth equality is the replacement of sequences of spatial diffeomorphisms (D\f n \...D\f2]D\fi\) by their product, 
which is (£%„]). The penultimate equality is the key one; it follows directly from equation (|l5|), namely from the 
transformation properties of the hamiltonian constraint under spatial diffeomorphisms. The last equality follows again 
from the definition of ordered exponential. In words, we have shown that the temporal evolution generated by the 
Lapse and the evolution generated by the Shift can be disentangled. 

While the manipulations above are formal (they are made inside a limit), the result itself is geometrically obvious: 
we can always rearrange the coordinates so that the Shift is zero, and compensate with a finite change of space 
coordinates at the end. If we do so, the Lapse N must be replaced by the Lapse in the new coordinates, which is N^. 

If the lapse is constant in space, = N. Then Un,o can be expanded as 



U N0 = 1 + (-») / dt C[N(t)] + (-if / dt / dt' C[N(t')] C[N{t)] + ... . (45) 
J o Jo Jt 

Its matrix elements between two spin network states can be expanded as 



(S f \U Nl0 (T)\Si) = (S f \Si) + (-») / dt(S f \C[N(t)]\Si) + (-if dt dt' (SflCWt')}^) (Si\C[N(t)]\Si) 

Jo Jo Jt 



(46) 



where we have inserted a complete set of intermediate states \S±) (S±\ (over which summation is understood). Using 
the explicit form (|l3|) of the hamiltonian constraint operator, we have 



(5/lf/jv.ol^) = (S f \Si) 

+H)J q dt [ J2 N(t, Xa )A a (Si) (Sf\D a \S t )+ N(t,xp)MSf) (Sf\Dl\Si 



-if I dt I dt' ^ ^ N(t,x a )N(t',x a ,) A a (Si)A a ,(Si) (S f \D a \Si)(Si\D a ,\Si) 

* a£[Si] a'e[Si] 







+ ... (47) 

(the second order term has three more summands, corresponding to the DD',D'D,D'D' terms). The first point to 
be noticed in this expression is that the sum over the intermediate state Si is finite. This is because both D and 
yield a finite number of terms only, when acting on a spin network state.^J Thus, the above expression is finite order 
by order. Next, the integrations can be performed explicitly, using (]34|). We obtain 



8 Lewandowski P8| has noticed that this finiteness might fail because of the moduli parameters of high valent intersections 
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(S f \U\Si) = (S f \Si) 



-HT) Yl A ^ (Sf\D a \Si}+ MSf) (Sf\Di\Si) 

\ae[Si] ae[S f ] 

-^P~ J2 E Aa( - S ^ Aa/( - S ^ {S f \Da,\Si) (Si\D a \Si) 

ae[Si] a/e[Si] 



(48) 



The structure of the expansion is now rather clear. At each order n, we have the D operator acting n times, n factors 
A, and a finite number of terms, coming from summing over nodes,hnks and e = ±1. 

Our next step is to integrate over Shift and Lapse (satisfying (pq)). The integration over lapse is trivial, as its 
dependence has dropped out the integral. This confirms the independence from the lapse that was mentioned in the 
previous section. The integration over the Shift amounts to imposing the Diff constraint. Indeed, it turns out to be 
equivalent to an integration over the diffeomorphism group, as in the group integration technique for solving the Diff 
constraint. Using (|44|) , we have 

U(T) = AT / [dN] ( [dN] D[g[N]] U Nfl , (T), (49) 



where we have explicitly indicated the dependence of g on N for clarity. We change integration variable N — > N^ 
(the Jacobian must be one, since this amounts to a change of coordinates )), and obtain 

U(T) = fif> (j [dN] U Nfi {T)\ (| [dN] D[g[N]]\ . (50) 

The TV integration can be traded for an integration over Diffo changing variables from TV to g[N], so we obtain 

U{T)=N ( [dg] D[g] U Nfi {T) (51) 

J Diffo 

for an arbitrary (irrelevant) choice of N satisfying (|34|), say N = T. The matrix elements of this operator are given 
by 



(S J \U{T)\S i )=M I [dg] (g.S f \U N , (T)\Si). (52) 

J Diffo 

The operator U(T) is now well defined in TLdiff^- Indeed, it is immediate to see that it is diff-invariant. For every 
two s-knots Si and Sf in Tidiff, we can arbitrarily pick Si and Sf such that Si € and Sf G Sf, and we have the key 
result that 

(s f \U(T)[si) = {S f \U(T)\S t } (53) 

is well defined (independent from the Si and Sf chosen). 

Furthermore, the operator D, depends on an arbitrary regularization -the location of the added link-, but a 
moment of reflection shows that the dependence on the regularization drops out in the step from U N fi(T) to U(T), 
by integrating the Shift. The reason is that different regularizations are related to each other by a finite diffeomorphism: 
the states D a \S) and D' a \S), where D and D' indicate two different regularizations of D are in the same s-knot: their 
difference becomes irrelevant in the scalar product (^3|). This result is due to the fact that all the factors in the 
expansion are individually well defined at the diffeomorphism invariant level. More precisely we have that 

(s f \ J2 Msi)D a \si) (54) 

a£[si] 



which were studied in |49|. The role of these parameters in the theory, however, is unclear. Finiteness of the proper time 
expansion may indicate that the correct version of the theory is the one in which the moduli parameters are removed, as 
suggested by many, and recently detailed in M . 
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is not only well defined, but also independent from the regularization of D. This fact allows us to write our expansion 
directly in diff-invariant form as 



(s f \U(T)\ Si ) = (s f \ Si ) 

H-(-iT) J2 A ^ (sf\D a \si}+ A «( s /) < s /l^h) 

\ae[si] ae[s f ] J 

+ 2! 51 51 Aa ( Si ) A a>(si)(sf\D a> \si)(si\D a \si) 

a£[si] a/e[si] 

+ ... . (55) 

This expression gives the 3-geometry to 3-geometry propagator of quantum general relativity as a series finite at 
every order. Notice that the expansion is in power of T, or, equivalently, in inverse powers of the Planck length, 
because this must divide T in order to recover physical dimensions. The utility of a perturbation expansion in inverse 
powers of G has been advocated in quantum gravity by Isham, Teitelboim and others plj ]. Intuitively, we can think 
that this quantity represents the probability amplitude that if we have a quantum state of the gravitational field (a 
quantum 3-geometry) \st) over a surface E^, we will find the quantum 3-geometry |s/) on the surface E/ in a proper 
time T . 



IV. SUM OVER SURFACES 



Surfaces in spacetime provide a natural bookkeeping device for the terms of the expansion (|5q ) in the same manner 
in which Feynman graphs provide a bookkeeping device for conventional QFT perturbation expansion. This fact leads 
us to give a nice graphical interpretation to the expansion ( |55| ) . 

Consider the 4-d manifold M. = [0, 1] x S. Denote the two connected components of the boundary of Ai as E^ and 
E f . We now associate a 2d colored surface a in M. -defined up to 4d diffeomorphisms- to each nonvanishing term of 
the sum in the right hand side of ( |55| ) . We begin by drawing s; in E; and s/ in E/ . As Si and s / have no information 
about the actual location of the graph, location is chosen arbitrarily (that is, up to a diffeomorphism in Diffo). The 
first term in (55) (zero-th order in T) is nonvanishing only if Sf = Si. In this case, let us slide Si across M. from E^ 
to E/, in such a way that it ends up over s/. To the term of order zero we associate the surface a = Sj x [0, 1] swept 
by Si. See Figure 2. 















\ / 

/ 

/ 











FIG. 2. Surface corresponding to a term of order zero. 

Notice that this is possible because the two spin networks are in the same s-knot. The surface we obtain is formed 
by 2d faces -submanifolds of M- joined along edges. The faces are swept by the spin network links, and the edges 
are swept by the spin network nodes. We color every face with the color of the corresponding link of Sj, and every 
edge with the color of the corresponding node of Sj. 

The surface associated to one of the summands of the second term (first order in T) in J55| ) is then defined as 
follows. In each summand, one of the nodes of Sj, say the node i, is altered by the operator D. st has two nodes 
more than Sj, say i' and i". We begin by sliding Si into the manifold by an arbitrary finite amount, until a position, 
say, s. Let p be the point in which the node i ends up. Then we slide s/ from E/ through the manifold in such a way 
that it converges to s. The three nodes and i" of s/ converge all three to p. We obtain a surface a, bounded by 
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Si and Sf formed by faces that meet along edges; four of these edges meet at the point p. We call p a vertex of the 
surface a. At the vertex p, a branches. Notice that four edges and six faces meet in p. See Figure 3. 





V 
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/ 

1 p / 

)/ 













FIG. 3. Surface corresponding to a first order term. 

We can imagine M. as a spacetime and Si as evolving continuously in a coordinate t from Si to Sf. At the spacetime 
event p, the spin network branches: the node i generates the two new nodes i' and i", which are born at i and then 
move away. A new face, spanned by the new edge that joins i' and i", is born in p. The branching represents the 
elementary vertex of the theory and is represented in Figure 4. 




FIG. 4. The elementary vertex. 



The generalization of this construction to higher terms is immediate. A term of order n in T corresponds to a 
surface a with n vertices. The (time) order in which the n D-operators act determines an ordering for the vertices. 
An example of a term of order two is given in Figure 5. It represents the transition from the s-knot with two trivalent 
nodes connected by three links colored (3, 5, 7), to the s-knot with the same graph, but colored (3, 6, 8). The 
intermediate step is the s-knot si, with four nodes. 
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FIG. 5. A term of second order. 



In the construction we have described, each vertex has 4 adjacent edges. Some of these edges are generated by 
the nodes of the incoming s-knot (the one at the right of the operator D corresponding to that vertex) and some by 
the outgoing s-knot. At each node, we denote the first ones as "past" edges and the second ones as "future" edges. 
Thus, each edge emerges as a future edge from one vertex, or from the initial hypersurface, and ends as a past edge 
in another vertex vertex, or in the final hypersurface. This defines a partial ordering of the vertices of each surface. 

A short reflection will convince the reader that all the surfaces that we obtain satisfy the following property. Each 
face has the topology of a disk, and the ordered set vertices around a single face has at most one local maximum and 
most one local minimum. We say that a colored surface a is "well-ordered" if all its faces satisfy this property. 

Now, observe the following, (i) The colored surface a (with the vertices ordered) contains all the information needed 
to reconstruct the corresponding term in the expansion (|55|). In fact, the factors A a (s) depend only on the coloring 
of the surface, (ii) Any well-ordered branched colored surface er, with colorings satisfying Clebsch-Gordan conditions 
at the edges can be obtained from a term in ([55]) . (iii) Two surfaces correspond to the same term if and only if there 
is a 4-d diffeomorphism that sends one into the other. 

These facts allow us to rewrite the expansion ( |55| ) as a sum over diffeomorphic inequivalent well-ordered surfaces a 
bounded by s,; and Sf. Therefore we can write the propagator (see (|35|)) as a sum of terms labeled by topologically 
inequivalent branched well-ordered colored surfaces^] s bounded by initial and final state: 

P(s f , Si ;T)= J2 ( 56 ) 

tr 

da— Si Us f 

The weight *4[<t](T) of the surface a is given by a product over the n{&) vertices of a: 

A[a](T)= [ - > JT A v (a). (57) 

The contribution A v (a) of each vertex is given by the coefficients of the Hamiltonian constraint defined in equation 



9 For a general description of such surfaces and their properties, see for instance [B2J and references therein. 
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A v (a) = A a (s). 



(58) 



(The nonvanishing matrix elements {s'\D a \s) have value one.) 

Can we attribute a physical meaning to the surfaces that enter the sum? The answer is yes. There is a natural 
way of interpreting a branched colored surface a as a discrete ("quantum") geometry. This geometrical interpretation 
was proposed in j| in a slightly different context; it holds in the worldsheet formulation of the simplicial model of 
GR 1 1] . First of all, consider a triangulation T of the manifold M. , and assume that the surface a sits over the dual 
two-skeleton of the triangulation. As we shall see in Appendix B, this is the natural way of viewing the surfaces a. 
Let a triangle (two-cell) S of the triangulation T be punctured by the face / (say with color p) of a in a point. Now, 
recall that according to canonical loop quantum gravity the colors of the spin networks are quanta of area: the area 
of a surface S pierced by a single link with color p = 2j is EG] 



A(S) = 16nhG VJ(J + rj. (59) 

In the spacetime picture, a link sweeps a 2d face /, which intersects S at a point. It is natural to suppose that the 
area of any spacetime 2-surface S is similarly determined by the coloring of the worldsheet. For instance, we may 
consider a three-dimensional hypersurface S that contains S, view the intersection between £ and the colored surface 
a as the "instantaneous position of the spin-network state on the ADM time X" , and assume that the results of the 
canonical theory can be applied. If we make this assumption, then we can say that the area of S is A(S) given in 
([59|). Therefore, a surface a assigns a (possibly vanishing) area to each triangle of the triangulation T . But fixing the 
areas of the triangles of a four dimensional triangulation is equivalent to fixing a discretized 4-geometry. Assigning 
areas is analogous to assigning the lengths of the links of the triangulation as in Regge calculus^ Thus, a surface 
a defines a discretized 4-geometry. The idea that areas of triangles could be variables more suitable than lengths 
of links in 4 dimensions was considered in ]l^ , |53| . Finally, more in general, we can say that a natural geometrical 
interpretation of the colors associated to the faces is the following: if a face has color 2j, it contributes a quantum of 
area 16whG \J j{j + 1) to the area of each spacetime 2-surface S at each point where it pierces S. 

This geometrical interpretation is "natural" , but not necessarily correct. In particular, the relation between the 
proper time T, and the spacetime geometry defined by the colors of the worldsheet is not clear. This relation should 
be investigated before taking the geometrical interpretation too seriously. 



A. Reconstruction of the A a (S) coefficient from surface data 

The coefficients A a (s) can be reconstructed directly from the colored surface as follows. Let a vertex v have n; past 
edges and rif future edges. A v (a) is non vanishing only if rii = 1 and rif — 3 or if rii = 3 and rif = 1. In this case, 
A a (s) is determined by the matrix elements of the hamiltonian constraint. 

It is very instructive to give an explicit construction of A v (a). Consider a 4d neighborhood B of the vertex v. 
Consider the 3d boundary dB of B. Let S v be the intersection between a and dB. A short reflection will convince 
the reader that S v is a colored graph in the 3-d space dB, having 4 = rii + rif nodes (that satisfy Clebsh Gordan 
relations), corresponding to the intersections between the 4 edges emerging from p and dB. See Figure 6. 





FIG. 6. The construction of the spin network Sv from the intersection of the surface with the boundary of a 4-sphere 
surrounding the vertex. The spin network is then cut into its past and future components S Vl i, S v j. 



There is a difference: in order to define a geometry, the lengths of the links must satisfy certain inequalities. The areas of 
the triangles must satisfy certain inequalities, as well as some equalities among them. Namely they are not all independent. 
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Now, cut all the links of S v that go from a past node to a future node. This procedure breaks S v into two spin 
networks with (equal) open ends, which we denote as S Vl i and S v j. (More precisely, these are s-knots, because 
they are determined only up to diffeomorphisms.) The value of the vertex is given by the matrix elements of the 
Hamiltonian constraint between these two spin networks. Namely, 

A v (a) = (S vJ \C[l]\S V!i ); (60) 

see (|ll|). This expression gives A v (a) as a function of the colorings of the edges and faces adjacent to the vertex 
p. This function is universal, and characterizes general relativity, in the same manner in which the Feynman vertex 
factor characterizes a QFT. We compute the vertex function A V (S) explicitly for the simplest case in Appendix B. It 
turns out to be expressed in terms of SU{2) n-j symbols of the colorings. 

Notice that it is the locality of the hamiltonian constraint that allows the sum over surfaces construction. This is 
a peculiar form of locality, quite different from conventional QFT locality: The action is not local with respect to a 
background structure, but with respect to the spin networks themselves. 

Equations (56-6^) provide a definition of the proper time propagator of quantum general relativity as a topological 



sum over branched colored surfaces. They represent our main result. 

V. CROSSING SYMMETRY 

Above, we have considered a reformulation of loop quantum gravity as of a sum over surfaces. Here we propose a 
modification of the theory, suggested by the reformulation. 

The value of the vertex A v {a) that we have computed in the last section depends on two inputs. First, on the 
coloring of the edges and faces adjacent to the vertex v. Second, on the distinction between "past" and "future" 
edges, namely on the way the vertex is located and oriented within the surface a. We suspect that the appearance 
of this orientation dependence is a sign that something has got wrong in the definition of the theory The action 
of GR is local and 4-d diff invariant, and therefore the action of the 4-geometry of a small region ("the vertex") is 
independent of how this region is sliced by equal time slices .0 

Thus, we propose a modification of the theory in which the orientation dependence is removed. We say, in general, 
that in a theory defined by a sum over branched colored surfaces, with weights given by products of vertex factors, 
the vertex is "crossing symmetry" if its value depends on the adjacent colorings only, and not on the distinction 
between past and future edges. BF theory |ll[] and simplicial GR H|J are theories of this kind, and have crossing 
symmetric vertices. In section V.A., we study the modification of the geometry of the vertex required to make it 
crossing symmetric. 

We then say that a hamiltonian constraint operator H has crossing symmetry if it defines a crossing-symmetric 
vertex via eq. (|6^) . The modification of the vertex that we consider in section V. A might be obtained from a different 
factor ordering of the hamiltonian constraint, and, as we shall see below, is strictly related 4d diff-invariance. Thus, 
here we are exploring the idea that 4d diff-invariance might fix residual factor ordering ambiguities. Of course, it 
should not be surprising that a spacetime formalism could simplify the discussion of 4d diff-invariance, a notoriously 
tricky issue in the hamiltonian framework. 

Let us make clear that we present crossing symmetry only as a proposal to be explored. We do not have a rigorous 
derivation of crossing symmetry from first principles, but only a heuristic plausibility argument, which we better 
detail below. 

Consider the path integral that formally defines U(T) in a proper time gauge, namely in a gauge in which the lapse 
is spatially constant. Consider a 4-metric g, that contributes to this path integral and a small spacetime region R, and 
let gn be the restriction of g to R. The region R is sliced by the proper time slicing. Let An(g) be the exponential of 
the action of this region. If the region is small enough, we can think of A^{g) as the matrix element of the evolution 
operator between "before i?" and "after R" , where "before" and "after" are determined by the proper time slicing, 
and thus identify An(g) with the vertex A v (cr). Now consider a different 4-metric g' in the integral, containing a 
region i?', such that g' R , is isometric to gn, but sliced in a different manner by the proper time slicing (the reader 
will easily convince himself that such a metric exists in general). Since the action is local and 4-d diff invariant, the 
contribution of gu to the sum must be equal to the contribution of g' RI , namely An(g) = ARi(g'). This implies that 
the matrix elements of the proper time hamiltonian between "before" and "after" according to one slicing of R ought 
to be the same as the matrix elements of between "before" and "after" according any other slicing. In other words, 



L We are dealing with the Euclidean theory, so there is no light-cone structure that defines local notions of past and future. 
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the matrix elements should be invariant under a 4d rotation of R that changes what is before and what is after. If 
we require the same to hold in our sum over surfaces, we obtain the requirement that vertices be crossing symmetric. 

This discussion shows that there is a relation between 4d diff invariance and crossing symmetry, because a 4d 
diffcomorphism "rotates" the vertex in 4d. Recall that the 4d diff invariance of the classical theory is expressed by 
the Poisson brackets 

{C[N} , C[M] } = C[NdM - MdN] . (61) 

One of the hard problems of the hamiltonian quantization program is to define a quantization of the hamiltonian 
constraint yielding a 4d diffeo invariant quantum theory. In particular, implementation of 4d diff invariance is 
presumably the missing ingredient for fixing quantization ambiguities of the hamiltonian constraint. Recall that the 
ambiguity in the definition of C[N] was fixed in |32) and J|3| to a large extent arbitrarily. Full implementation of 
the quantum version of (|l]) should ensure 4d diff invariance, but has proven hard to realize. We are therefore lead 
to the suggestion that we can cure the slicing dependence by taking advantage of the remaining operator ordering 
ambiguity, and at the same time, cure the excessive ordering ambiguity by imposing slicing independence. In other 
words, we can impose some form of J^d diff -invariance requirement in order to reduce quantization ambiguity. Here we 
are suggesting that in a covariant formalism crossing symmetry might be the key for implementing 4d diff invariance. 

A. Vertices with crossing symmetry 

Let us investigate meaning and consequences of requiring crossing symmetry. First, we should require that different 
ways of cutting S v with three nodes on one side and one node on the other yield the same A v (o~). This leads to 
conditions on the A a (S) coefficients, that will be studied elsewhere. A more interesting case is the following. First, 
let us help intuition by redrawing the elementary vertex of the theory (Figure 4) in a more symmetric way. This is 
done in Figure 7. (For simplicity, we restrict the following analysis to trivalent nodes.) 




FIG. 7. The elementary vertex and its associated spin network S v . 



There are five topologically inequivalent ways of cutting S v , giving, respectively (rii, nj) (number of initial and final 
nodes) equal to (0,4), (1,3), (2,2), (3,1) and (4,0). The last two are the time reversal of the first two, leaving three 
genuinely independent cases. In Figure 8, we show the possible cuts, and the corresponding spin networks transitions 
in the hamiltonian picture. Time reversed cuts give just the opposite transitions. 
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FIG. 8. The (1, 3), (0, 4) and (2, 2) cuts of the elementary vertex, and, in the last column, the corresponding spin network 
transitions in the hamiltonian picture. (For (3,1) and (4,0), look at (1,3) and (0,4) upside down, and reverse the arrow of the 
transition.) 

Case (1,3) is the one described in the previous section (Figure 6). Crossing symmetry requires that the hamiltonian 
constraint generate the transitions (0,4) and (2,2) - described in the last column of Figure 8- as well, with the same 
amplitude. 

Consider these two new transitions. We begin with (0,4) (second line in Figure 8). This transition represents a 
matrix element of a hamiltonian that creates a "small" tetrahedron from the state with no loops. The fact that 4d 
invariance requires the presence of such "birth" terms has already been argued, on general grounds, in |54| . In terms 
of surfaces, the term looks as in Figure 9. 




FIG. 9. The (0,4) transition: creation of a tetrahedron. 

This is the very same surface as in Figure 4, o,nd in Figure 7; but drawn with a different orientation in "spacetime" . 
With this orientation, it describes a tetrahedral spin network emerging from nothing. 

Can such a term originate from an ordering of the hamiltonian constraint? Surprisingly, the answer is positive. 
We sketch here an hand waving argument. The regularized hamiltonian constraint (FEE) is formed two parts: s 
"small loop" that corresponds to the classical curvature (F) term; and the term that "grasps" : the two hands of the 
T 2 operator in |lj],|32| (or the Volume operator in Q), corresponding to the triads EE (or the triads multiplied by 
a suitable density factor). Traditionally, the order chosen is FEE: the "small loops" is added after the grasping. 
Reverse this order, choosing EEF, and have the small loop being inserted first. Then the hamiltonian constraint has 
non- vanishing action on the vacuum as well, because the grasping term can grasp the "small loop". In particular, 
this may create a "small" tetrahedron. For instance, in the construction in 
precisely, a tetrahedron. See Figure 10. 




FIG. 10. Creation of a tetrahedron from the self grasping of T 2 . 



17 32], the T 2 can grasp itself, producing, 
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The (2,2) term (third line of Figure 8) gives a rearranging of two nodes. The corresponding surface looks as in 
Figure 11. 




FIG. 11. The (2,2) vertex. 
Again, this is just a different orientation of the same elementary vertex. 

Now, we could search for an ordering of C[N] yielding a vertex having crossing symmetry. But this task is 
superfluous since we already know what we should obtain. We can directly postulate that the hamiltonian constraint 
yield crossing symmetry, and deduce the amplitudes of the (0,4) and (2,2) matrix elements from the value of the (1,3) 
vertex.^ 

There is a physical motivation that supports the above argument. It has been observed ]5f| that the orderings of 
the hamiltonian constraint studied so far generate a dynamical evolution that appear to be excessively "local" . They 
preserve the general structure of the network on which they act, simply "dressing" nodes. These difficulties have been 
recently detailed in Ref. ||, where it is argued that no long range interaction is likely to emerge from a hamiltonian 
with these features. As pointed out by Thiemann, the argument is far from conclusive, because it contains a jump 
from the non-physical (gauge) coordinate evolution to the physical one, and this jump may be un ungranted. In any 
case, adding the new vertices (2,2) and (0,4) would cure these potential difficulties. 

Also, we note that one of the consequences of adding the new vertices is that the faces of the surfaces in the sum 
do not need anymore to be (topologically) two-disks, as followed from the original Hamiltonian expansion, and there 
is no sense in the requirement of the surfaces being well-ordered. 

Finally, notice that one could also search for the form of the vertices from general a priori requirements. An arbitrary 
crossing symmetric vertex is obtained by replacing ( |6fj| ) with a function A(S V ) of the spin network S v associated to 
the vertex 

A v (a) = A(S V ). (62) 

For instance, for trivalent nodes, S v is a tetrahedron: the function A( ) must respect tetrahedral invariance. Notice 
that there are not many functions with these features. A natural choice is 



A(S V ) = Tet 



a b c 
d e / 



(63) 



where a — / are the colors of the links of S v and Tet is the totally symmetric form of the 6-j symbols ( |36| and 
p3|). We think that a theory defined in this way is worth exploring. 



VI. CONCLUSION 



Our main result is contained in the equation 



P(s f , Si ;T)= £ ( 64 ) 



ve\cr\ 



f 



12 Notice that the symmetrization of the hamiltonian constraint given in (g|) can be seen as a first step in this covariantisation 
of the operator: it is equivalent to the postulate that the (3,1) cut has the same value as the (1,3) cut. 
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[see ( p6[p7| )], which expresses the dynamics of quantum general relativity in terms of a sum over surfaces a. 

More precisely, the proper time propagator of quantum GR can be expressed in terms of a sum over topologically 
inequivalent branched colored surfaces, bounded by the initial and final s-knots. The contribution of each surface to 
the sum is the product of one factor per each vertex (branching point) of the surface. The contribution a each vertex 
is a simple S?7(2)-invariant function A v (a) of the colors of the faces and edges adjacent to the vertex. This function 
characterizes the quantum theory in the same manner in which the Feynman graph vertices characterize a quantum 



field theory. The vertex A v (a) of general relativity is given by a product of Wigner in — j symbols 1 34 . 

The essential property of the expansion (|64|) is that it is finite order by order, and explicitly computable. This 
finiteness is intriguing. In order to calculate physical quantities, we must have the proper time propagator for 
multifingered proper times and we must integrate over the multifingered proper time. We expect that the integration 
could yield finite results if performed over expectation values of appropriate physical quantities. Work is in progress 
in this direction, and will be reported elsewhere. 

We close with some comments. 

• Our construction is strongly reminiscent of discretized quantum gravity on a lattice ||,|^,|8U60), particularly 
in its simplicial formulations ||,[53| . It is shown in Q that one can discretize general relativity over a simplicial 
lattice, and express the gravitational degrees of freedom as colored branched surfaces over the (dual) two- 
skeleton.^] Even more remarkably, the partition function is given in the discretized case by a construction very 
similar to that given here: the contribution of a vertex is determined by the intersection between the boundary 
of a 4-simplex around the branching point of the surfaces, and the surface. This defines a spin network S v , 
which, in the discretized case, can be any subgraph of the 1-skeleton of a 4-simplex. Therefore vertices have up 
to 5 edges and 10 faces (see Appendix B) in the discretized case. In this paper, nonvanishing vertices have 4 
edges and an arbitrary number of faces. Thus, the simplicial construction corresponds to a cut of the sum ( |56| ) 
in two respects: the maximum number of vertices is fixed by the triangulation, and vertices have 10 faces at 
most. 

• One can view the sum over surfaces defined here as a version of Hawkings' integral over 4-geometries. Indeed, 
a colored two surface defines a discrete 4-geometry. The integral is replaced here by a sum, and explicit 
computation can be performed. Presumably, the construction then be used to define a number of related 
theoretical tools as partition function, Hartle-Hawking state, and similar Q . 

• Each individual term in the expansion (^5|) is finite. Divergences can arise in summing the series, and in 
integrating over proper time. 

• The similarity with the formulation of string theory as a path integral over worldsheets is tantalizing. On this, 
see the discussion in [Q. The dynamics is different. In string theory, the contribution of each surface to the sum 
is given by the area of the surface, and therefore it depends on a fixed background metric on the manifold. Here, 
on the contrary, the contribution of each surface depends only on the (coloring and) topology of the surface. 
Thus, quantum GR resembles a "background independent" version of string theory. The techniques developed 
here could perhaps have relevance for connecting loop quantum gravity with string theory ]7|j62]| - or for the 
construction of a non-perturbative background independent formulation of string theory. 
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1 Surfaces seem to be playing an increasing role as a way to capture the gravitational field degrees of freedom. See for instance 
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APPENDIX A: TERMINOLOGY 



To help the reader, we collect here a list of terms employed. 
Node: Point in 3d space where the links of a spin network meet. 
Link: Line in 3d space connecting two nodes of a spin network. 
Face: Surface in 4d spacetime (swept by a link). 

Edge: Line in 4d spacetime where several faces meet (swept by a node). 

Vertex: Point in 4d spacetime where several edges meet. 

A spin network is formed by nodes and links. A branched surface is formed by faces, edges, and vertices. For 
the branched surfaces that live on the 2-skeleton of the dual triangulation of the manifold in simplicial BF theory, 
faces, edges, and vertices live on 2-, 1- and 0-cells respectively of the cellular decomposition dual to the simplicial 
triangulation. They are therefore associated to 4-, 3- and 2-simplices of the triangulation, respectively. Therefore a 
vertex corresponds to a 4-simplex, an edge to a tetrahedron, and a face to (its dual) triangle. 



APPENDIX B: COMPARISON WITH THE OOGURI-CRANE-YETTER 4D TQFT 

The structure of quantum general relativity in the form presented in this paper is surprisingly similar to the Ooguri- 
Crane-Yetter (OCY) four dimensional topological quantum field theory Jl(],[ll| a rigorously defined simplicial lattice 
version of four dimensional SU (2) BF theory. 

More specifically, our expression for the proper time propagator U (T) of GR as a sum over worldsheets resembles 
in many ways the worldsheet sum |§|5(|0 for the projector on physical states, Uocy, of the OCY model. 

The OCY model is a 4d generalization of the Ponzano-Regge-Turaev-Viro (PRTV) model ||l2] , |l3|| , which, in turn, 
can be seen as a quantization of 3d GR, or a quantization of 3d Chern-Simon theory. In it was shown that the 
PRTV model is a theory of the dynamics of spin networks (loops in the terminology of having the same physical 
interpretation as the spin network basis states in continuum 3d GR. Thus one might expect a similarity between the 
kinematic of the OCY model and loop quantized GR JlJ]. On the other hand, the 4d OCY model like the 3d PRTV 
model, but unlike 4d GR has no local degrees of freedom, so one also expects large differences between the theories. 

In this section, we sketch the OCY theory, outline a construction of the worldsheet sum for the partition function 
ZbF of the OCY model along the lines of §J^| and discuss its similarities and differences with the formulation of 
quantum GR presented here. We believe that this comparison helps illuminate the much debated issue of the relation 
between quantum gravity and TQFTs Jl5| ]. 

We introduce here Ooguri's original version Jh} the OCY model heuristically, as a discretization of BF theory 
without cosmological constant. BF theory is given in terms of two fields, an SU(2) connection A 1 , with curvature F % , 
and an sit (2)-algebra valued 2-form B\ by the action [B3| 



Sbf = J B l AF\ (Bl) 

Before proceeding, it is interesting to note that conventional general relativity can be obtained from BF theory by 
simply adding a constraint term. Indeed one can show that the theory 

Sgr, = I B l AF l + cf> ij B i A B 3 , (B2) 



where the Lagrange multiplier <f>ij is traceless and symmetric, is equivalent to general relativity 64 1. 
Consider the partition function of the BF theory 



14 The worldsheet sum of [|56j is actually for the 3d Ponzano-Regge-Turaev-Viro (PRTV) model, but is easily extended to the 
OCY model. Iwasaki M] has proposed an interesting alternative, a closely related worldsheet formulation of the PRTV model 
which is also easily extended to the OCY model. See also 

15 The construction outlined here is a sort of baby version of that for simplicial GR in |54|]. 



21 



Zbf = J [dA] [dB] e _i / fl * Ajrt . (B3) 

Integrating over B, we obtain 

Z BF = J[dA] 6[F], (B4) 

namely an integral over flat SU(2) connections. Let us define a lattice version of this theory by fixing a simplicial 
decomposition of the 4d manifold. (See also f(5(J.) 

Consider the dual of the simplicial decomposition. There is one element of this dual cellular decomposition that 
plays a central role in the construction: the "wedge". Consider a dual-two-face /. A dual-2-face is a 2d polygon. It 
intersects a two-face of the simplicial decomposition in a "central" point o. Its vertices are centers of 4-simplices and 
its sides are lines connecting such centers. Each of this sides, which connects the centers of two simplices, crosses 
the tetrahedron that forms the boundary between the two simplices. Let p be the crossing point. Each point p cuts 
one of the sides of the polygon /. By drawing lines connecting the the points p to O, we divide the polygon / in 
quadrangles, called "wedges" . A wedge is thus a 2d quadrangle that has four sides: two of these are Id lines that join 
centers of 4-simplices with (the center p of) a bounding tetrahedron; these are denoted 3-4 flags. The other two lines 
join the center of a tetrahedron with the center of a two-face. These are denoted 3-2 flags. See Figure 8. 



4-symplex 




Wedge 

FIG. 12. A dual-two-face and its decomposition in wedges. 



We choose to represent the connection by means of group elements associated to Id elements in the dual cellular 
decomposition. More precisely, we associate a group element U to each 2-3 flag (segment connecting the center of a 
2-simplex with the center of one of the tetrahedra surrounding it), and one group element W to each 3-4 flag (segment 
connecting the center of a tetrahedron with the center of an adjacent 4-simplex). These group elements can be thought 
as the exponential of the connection along the segments. Each wedge w is bounded by four of such segments (two 
of the 2-3 kind and two of the 3-4 kind); let Ui(w), 1/2(10), Wi(w), W2(w) be the group elements associated to the 
segments that bound the wedge w. We can express the requirement that the conne ctio n is flat by requiring that the 



holonomy of the connection around each wedge is trivial. Then a discretization of (B4) is given by 



Zbf = f[dU] [dW] Y[S[U 1 (w)U 2 (w)W 1 (w)W 2 (w)} (B5) 
J in 

where the delta function is the delta function of the unit on the SU (2) group. We can expand the delta function in 
characters. For each w, we have 

8[U l (w)U 2 (w)W l (w)W 2 (w)] = J2( 2 J + l)Tr j (U 1 (w)U 2 (w)W 1 (w)W 2 (w)), (B6) 

3 

where j labels the irreducible representations of SU(2), and Trj(U) is the trace of the group element U in the 
representation j. Using this, we can rewrite ( |B5| ) as an integral over group elements U and W associated to segments 
and half integers j associated to wedges: 

Zbf = f[dU] [dW]^2Yl(2j(w) + l)Tr j{w) (U 1 (w)U 2 (w)W 1 (w)W 2 (w)). (B7) 

We can view the group elements U and W as the discrete version of the connection A, the j's as a discrete version 
of the two form B and the expression Trj( w \(Ui(w)U 2 (w)Wi(w)W 2 (w)) as the discrete version of the expression 
exp(-iB l AF l ). 
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Next, let us perform the group integration in (B7) explicitly. By integrating over the group elements U we force 
the colors of the wedges belonging to the same dual 2-cell to be equal. By integrating over the group elements W we 
force the j's of the (four) dual 2-cells that join on a 1-cell to satisfy Clebsh-Gordan relation where they meet, leaving 
an extra degree of freedom J associated to each such dual 1-cell; J runs over the independent couplings of four SU (2) 
representations. Finally, we end up with numerical factors associated to the 0-cells of the dual triangulation (plus 
other factors associated to faces and edges, which we disregard here in order not to make the exposition to heavy). 
Such numerical factors turn out to be 15- j symbols associated to the five J's of the five 1-cells and the ten j's of the 
ten 2-cells adjacent to each vertex.^. Performing these integrations explicitly is a simple and interesting exercise. 
After these integration over the group elements, the theory is therefore reduced to a sum over colorings on the 2-cells 
and 1-cells, satisfying Clebsh-Gordan relations. We can interpret a zero color as no surface at all, and identify the 
2-cells with faces and the 1-cells with edges of branched colored surfaces. Thus, we can write the partition function 
as a sum over branched colored surfaces living on the dual two skeleton of the simplicial triangulation. We obtain 

Z BF = yA BF [a] (B8) 



where the contribution of each surface is (up to the face and edge factors we have disregarded for simplicity) a product 
of vertices's factors 

A B F[(r}=l[A BFtV {a); (B9) 

V 

the vertex factor is the 15-j symbol of the colorings adjacent to the vertex. 

The similarity of this result with the construction in this paper, equations ( |5^j57| , p8| ) , is striking. In both cases, we 
have a sum over the same kind of branched colored surfaces, and the weight for each surface is the product of vertex 
factors, where vertex factors are simple SU (2) invariant functions of the adjacent colorings. Thus, the structure of 
quantum general relativity turns out to be extremely similar to the structure of a topological quantum field theory. 
Of course there are differences, and these differences are crucial. Let us examine them in detail. 

• First of all, the worldsheet amplitudes in the BF theory that we are considering are the amplitudes in the 
projector on physical states, while the GR worldsheet amplitudes of the present paper are from the sum for 
the proper time propagator U(T), so we might be comparing apples and oranges. However, if we accept the 
not unreasonable hypothesis that U(T) is a partial sum of terms in a sum over surfaces for U in GR, we can 
compare the theories in a direct way. 

• The vertex factor is different in the two theories: in BF theory it is a 15-j symbol, while in GR it is a combination 
of 9-j and 6-j symbols. This difference depends on the different dynamics of the two theories, and should be at 
the root of the other differences. 

• In the case of BF theory there is a crucial theorem holding: triangulation independence. Refining the triangula- 
tion does not change the overall sum. This is the reason for which the theory is topological, and is a consequence 
of the fact that the classical theory has no local degrees of freedom. In GR, nothing similar holds, because GR 
has genuine local (although non-localized!) degrees of freedom. Therefore there is no reason to expect anything 
like triangulation independence for GR. 

• The ensemble of surfaces over which the sum is defined is different in the two cases. In the BF case, we sum 
over surfaces over a fixed triangulation. In the case of GR, we sum over all topologically inequivalent surfaces, 
with an arbitrary number of vertices. Therefore, in the case of BF theory the surfaces to be consider are finite 
in number. In the case of GR we have to sum over arbitrarily complicated surfaces, or, equivalently, sum 
over arbitrarily fine triangulations of the manifold as well. Notice that this difference is a consequence of the 
previous point, namely triangulation independence of BF. We could average over arbitrarily fine triangulations 
in BF as well; but this would not affect the result, because each triangulation yields the same contribution as 
the coarsest one. Therefore: Diff invariance of the sum is implemented in two different ways in BF and in 
GR, corresponding to the fact that BF is topological, while GR is not: in BF, invariance is obtained thanks to 
triangulation independence; in GR invariance is obtained by summing over arbitrarily fine triangulations. 



16 A dual 0-cell is always adjacent to five 1-cells and ten 2-cells, because it corresponds to a four-simplex of the original 
simplicial decomposition, which is bounded by five tetrahedra and ten faces. 
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• The rigorous version of BF theory requires SU{2) to be replaced by quantum SU(2). This can be seen simply 
as a smart stratagem for regularizing the sum in an invariant manner, yielding a finite result. Notice that 
in GR regularizing SU(2) to quantum SU(2) would not guarantee an overall finite sum, because the surfaces 
themselves are infinite in number. Thus, quantum GR does not admit a rigorous finite version as quantum BF, 
at least as far as we can presently see, even if one attempts to replace SU{2) with quantum SU(2) in GR [p5|. 

• The vertices of BF have always five edges and ten faces, while vertices of GR have (at least with the ordering 
considered so far) four edges and an arbitrary number of faces. 

These points illuminate the difference between quantum GR and topological field theories. Let us discuss this point 
in more detail. 

Both theories are invariant under diffeomorphism. However, diffeomorphism invariance does not imply that a 
quantum theory is topological in the sense of having a finite number of degrees of freedom. We expect GR to have 
an infinite number of degrees of freedom. Thus Atiyah's axioms for topological quantum field theory are likely to be 
suitable for quantum general relativity as well, if we drop the request that the Hilbcrt spaces attached to boundaries 
of the four-manifold be finite dimensional. 

Finally, we may turn the comparison the other way around, and describe quantum BF theory in the language so 
far used for quantum GR. We can capture quantum BF theory in terms of its vertex. A BF vertex has five edges and 
ten faces. Assume that one of these edges comes from the past, and four go to the future (the other cases are given 
by crossing symmetry, that clearly holds in BF theory) . A moment of reflection shows that the elementary vertex of 
BF theory "opens up" a four-valent intersection of a spin-network into a "small" tetrahedron: 




FIG. 13. The elementary vertex of BF theory. Notice that there are 5 nodes and 10 links, yielding 15 colors. The value of 
this vertex is the 15- j symbol of their 15 colors. 

The matrix element of the hamiltonian between these two (partial) spin networks is the symbol of the 15 colors 
associated to the four links and the one node of the incoming spin network, and the six other links and four other 
nodes of the newly created tetrahedron. It would be interesting to derive this hamiltonian from a hamiltonian loop 
quantization of BF theory. 



APPENDIX C: DIFF INVARIANT SCALAR PRODUCT 



We work out here an example of diff- invariant scalar product between s-knot states 40 1. (On the inner product 
between spin networks see [^7,Q, and |^3[ which we follow here). Let s be the s-knot defined as follows, s has three 
4-valent nodes, i, j and k, and the six links 



(ki,2), (kj,2), (fci,4), (fcj,4), (ij,3), (y,5), 



(CI) 



where each link is indicated by the two nodes it connects and its color. Explicitly: 




(C2) 



To specify the state, we have to give also the coloring of the nodes. We choose an expansion of the nodes in a trivalent 
graph (in each node) by pairing the two links colored 2 and 4. We have one virtual link for every node (we denote 
them ei, e,j and efc), which we assume to be colored as follows: ej : c, ej : 6 and : 2, where c (which can take the 
values 6, 4 and 2) will be specified later on 
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Next, let us define the s-knot s'. Let it be the same as above, but with a different coloring of the node j. We expand 
j by pairing (kj, 2) with and (kj, 4) with (ij, 3). Let the internal link have color 3: 

s' = i:W kW j\A (C4) 




Let us compute the scalar product (s|s'}. First, we have to list the automorphisms of the spin network (taking links' 
colors, but not node's colors into account). There is only one non trivial automorphism a: it exchanges i and j. Thus, 
(|) gives 

(s\s') = (S\S') + (aS\S') (C5) 
where S £ s and <S" £ s' have been selected to have the same graph, with the same colored links. The only contribution 



to (C5) comes from the nodes. The node k is the same in the two states and therefore gives no contribution (recall 
we have chosen normalized states). Thus, we have 

(s\s') = (S\S') t (S\S') :j + (aS\S') l (aS\S') :j (C6) 

where we have indicated with ( | ), the scalar product restricted to the space of one node. Spin-network states with 
the same trivalent expansion are orthonormal. The change of basis is given p3j by the recoupling theorem 

b \y c 

y< c =E{:2.;} a < c7 > 

ad' ad 
where the quantities | ^ ^ \ \ are su(2) six-j symbols (normalized as in (56)). This gives us immediately 



(s\s% 




(S\S')j 


: h 


[aS\S')i 


= s c6 


aS\S') 3 





5 3 6 



5 3c 



(C8) 
(C9) 
(CIO) 
(Cll) 



Therefore 



Thus, if c = 6 there are two contributions to the scalar product, one from each of the two elements of the automorphism 
group of the spin network and we have 

While if c = 2 or c = 4 

4 2 3 ] 56 



5 3 6 f " 75- ( C14 ) 
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